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ABSTRACT 



We develop and discuss global accretion solutions for viscous ADAF disks 
containing centrifugally supported isothermal shock waves. The fact that such 
shocks can exist at all in ADAF disks is a new result. Interestingly, we find 
that isothermal shocks can form even when the level of viscous dissipation is 
relatively high. In order to better understand this phenomenon, we explore all 
possible combinations of the fundamental flow parameters, such as specific energy, 
specific angular momentum, and viscosity, to obtain the complete family of global 
solutions. This procedure allows us to identify the region of the parameter space 
where isothermal shocks can exist in viscous ADAF disks. The allowed region is 
maximized in the inviscid case, and it shrinks as the level of viscous dissipation 
increases. Adopting the canonical value 7 = 1.5 for the ratio of specific heats, 
we find that the shock region disappears completely when the Shakura-Sunyaev 
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viscosity parameter a exceeds the critical value ~ 0.27. This establishes for the 
first time that steady ADAF disks containing shocks can exist even for relatively 
high levels of viscous dissipation. If an isothermal shock is present in the disk, it 
would have important implications for the acceleration of energetic particles that 
can escape to power the relativistic jets commonly observed around underfed, 
radio-loud black holes. In two specific applications, we confirm that the kinetic 
luminosity lost from the disk at the isothermal shock location is sufficient to 
power the observed relativistic outflows in M87 and Sgr A*. 

Subject headings: accretion, accretion disks - - hydrodynamics - - black hole 
physics — galaxies: jets. 

1. INTRODUCTION 

The accretion of matter onto compact objects is the fundamental mechanism power- 
ing a variety of high-energy astrophysical sources, such as low-mass X-ray binaries, massive 
black holes, and active galactic nuclei. Shakura & Sunyaev (1973) and Novikov & Thorne 

(1973) laid the foundation for our understanding of these sources by developing the first 
physical models for geometrically thin accretion disks. These early models did not treat 
the pressure and advection terms in the conservation equations correctly, and furthermore 
no attempt was made to satisfy the inner boundary conditions at the event horizon. In- 
stead, the disk was simply terminated at the marginally stable orbit. The inner regions 
of these disks were shown to be viscously and thermally unstable by Lightman & Eardley 

(1974) , and the model was subsequently improved by Paczyhski & Bisnovatyi-Kogan (1981) 
and Muchotrzeb & Paczyhski (1982) who incorporated advection and pressure effects into 
their models. Following a similar approach, the general global solution for accretion in the 
hydrodynamic limit, including advection, viscosity, and thermal effects, was obtained by 
Chakrabarti (1990, 1996). 

1.1. Advection-Dominated Accretion 

The ADAF model remains a central paradigm in contemporary black-hole accretion 
theory as a possible explanation for the dynamical structure of sources with significantly 
sub-Eddington accretion rates (see Yuan 2007 for a review). The gas in the inner regions 
of such disks is quite tenuous, and therefore the plasma radiates very inefficiently, leading 
to high temperatures. At large radii, the hot ADAF inner region is expected to merge with 
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a cool, geometrically thin Shakura & Sunyaev outer region characterized by a Keplerian 
angular momentum profile. Since any global accretion solution onto a black hole must cross 
the event horizon supersonically in order to satisfy the fundamental boundary conditions 
imposed by general relativity, ADAF accretion disks are necessarily transonic. Narayan, 
Kato, & Honma (1997) studied the behavior of transonic ADAF disks subject to the Shakura- 
Sunyaev viscosity prescription. The associated asymptotic boundary conditions governing 
the dynamics close to the event horizon, and the related physical restrictions on the viscosity 
prescription, were examined by Becker & Subramanian (2005) and Becker & Le (2003). 
Interest in the ADAF model remains strong, and it has been the focus of many recent studies. 
For example, Narayan & McClintock (2008) have analyzed the properties of ADAF disks 
close to the event horizon; Takahashi (2007) has explored the influence of causal viscosity 
prescriptions on the structure of ADAF disks; Ma, Yuan, & Wang (2007) have probed the 
role of magnetically-induced torques in ADAF disks; and Yuan, Ma, & Narayan (2008) have 
examined the properties of a simplified global model that utilizes an approximate form for 
the radial momentum equation. 



1.2. Critical Structure of Transonic Disks 

Due to the supersonic nature of the inflow at the event horizon, any physically achievable 
flow configuration must contain at least one critical point, where the flow transitions from 
subsonic to supersonic. However, for a given set of incident flow parameters (e.g., specific 
energy, specific angular momentum and viscosity), the flow may possess multiple critical 
points (Chakrabarti 1989a, 1989b, 1990, 1996; Chakrabarti & Das 2004; Le & Becker 2005). 
This suggests the possible existence of global flow solutions with a standing shock wave 
located between two critical points. In their early study of viscous, transonic ADAF disks, 
Narayan et al. (1997) focused on global solutions that pass through a single critical point, and 
consequently they did not obtain any shocked disk solutions. The problem of the structure of 
ADAF disks governed by the Shakura-Sunyaev viscosity prescription was recently reexamined 
by Becker, Das, & Le (2008), who demonstrated that isothermal shocks can exist in such 
disks even in the presence of substantial viscosity. 

In the present paper, we extend the approach taken by Becker, Das, & Le (2008) by 
performing a complete analysis of the structure of viscous ADAF disks containing either sin- 
gle or multiple critical points. Solutions with multiple critical points may include a standing 
shock, if it is possible to satisfy the shock jump conditions at some radius r* between the 
outer critical point (at radius r = r° ut ) and the inner critical point (at radius r = r™). 
Shocks occur when the accretion flow is impeded due to the presence of a "centrifugal bar- 
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rier" situated close to the horizon. Although such barriers can exist even in viscous disks, 
they begin to dissipate as the level of viscosity increases because of the increasing efficiency 
of the outward diffusion of angular momentum. Hence disks with large viscosity parameters 
are not expected to possess shocks, and our results confirm this. In the post-shock region, 
the flow accelerates towards the black hole, eventually passing through the inner critical 
point and crossing the horizon supersonically. 

1.3. Shock Waves in Viscous Disks 

The question of the stability of viscous accretion flows containing standing shock waves is 
currently not resolved. However, an important new development in the field was provided by 
Nagakura & Yamada (2008), who used both linear stability analysis and general relativistic 
hydrodynamical simulations to investigate the stability properties of shock waves in inviscid, 
advection-dominated accretion flows around Schwarzschild black holes. When subjected 
to strong perturbations around the initial condition, the authors find that the dynamical 
configuration is stable, in the sense that the shock continues to exist, although the shock 
radius may oscillate if the perturbations are non-axisymmetric. Similar results were obtained 
by Okuda et al. (2007). Although these authors do not consider the effect of viscosity in 
their work, their conclusions nonetheless support the general idea that shocks may be able 
to exist in a stable configuration in relativistic disks. Consequently their findings provide a 
sound theoretical basis for the study of the structure of shocked, viscous disks developed in 
the present paper. When a shocked solution is dynamically possible, we argue based on the 
second law of thermodynamics that a shock will form since shocked-disk solutions possess 
higher entropy than smooth solutions (Becker & Kazanas 2001). 

The study of shock waves in accretion flows around compact objects has been undertaken 
by a variety of authors using either analytical methods (Fukue 1987; Chakrabarti 1989a, 
1989b, 1990, 1996; Abramowicz & Chakrabarti 1990; Nobuta & Hanawa 1994; Yang & 
Kafatos 1995; Yuan, Dong, & Lu 1996; Caditz & Tsuruta 1998; Kovalenko & Lukin 1999; 
Das, Chattopadhyay, & Chakrabarti 2001; Das & Chakrabarti 2004; Le & Becker 2004, 
2005, 2007; Yu, Lou, Bian, & Wu 2006; Becker, Das, & Le 2008) or numerical simulation 
(Hawley et al. 1984a, 1984b; Chakrabarti & Molteni 1993; Molteni, Sponholz, & Chakrabarti 
1996; Ryu, Chakrabarti & Molteni 1997; Chakrabarti, Acharyya, & Molteni 2004). The 
physical configuration of the shock depends on the detailed microphysics of the dissipation 
mechanism operating at the shock. The three types of shocks most frequently studied are 
(1) Rankine-Hugoniot shocks, (2) isentropic shocks, and (3) isothermal shocks (Chakrabarti 
1990; Abramowicz & Chakrabarti 1990). Rankine-Hugoniot shocks possess the largest post- 
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shock temperatures because the energy flux is strictly conserved as the gas crosses the shock. 
Hence shocks of this kind are radiatively very inefficient. In an isentropic shock, the entropy 
generated via dissipation is immediately radiated away, and consequently such shocks display 
smaller temperature increases than the Rankine-Hugoniot variety. In an isothermal shock, 
the temperature does not increase at all because all of the dissipated energy and entropy 
escapes from the flow at the shock location. Hence isothermal shocks represent the most 
efficient mechanism for removing energy from the flow. 

Le & Becker (2004, 2005, 2007) demonstrated that the energy removed from the flow 
at the isothermal shock location can be understood physically as a consequence of the first- 
order Fermi acceleration of relativistic particles in the disk, combined with the diffusion 
and escape of the accelerated particles through the upper and lower surfaces of the disk. 
This approach facilitates the development of a completely self-consistent, coupled model 
that simultaneously describes the dynamical structure of the disk/shock system and also the 
energy spectrum of the relativistic particles escaping from the disk to form the observed jet 
outflows. We therefore focus on isothermal shocks in the present paper. 

The properties of isothermal shock waves in accretion disks have been studied previously 
by Abramowicz & Chakrabarti (1990), Lu & Yuan (1998), Fukumura & Tsuruta (2004), and 
Das, Pendharkar, & Mitra (2003). However, the studies cited above all focused on inviscid 
flow models, which probably do not adequately describe real accretion disks. In the present 
paper, we obtain for the first time detailed dynamical solutions for the structure of ADAF 
disks containing isothermal shocks. The remainder of the paper is organized as follows. 
In § 2 we discuss the ADAF model assumptions, the governing equations, and the critical 
conditions for the structure of ADAF disks. We discuss the isothermal shock jump conditions 
and examine global accretion solutions in § 3, and we explore the properties of the standing 
shock and classify the parameter space for shocked accretion solutions in § 4. Our dynamical 
formalism is used to model the development of the relativistic outflows in M87 and Sgr A*in 
§ 5, and concluding remarks are presented in § 6. 



2. GOVERNING EQUATIONS 

We consider a thin, axisymmetric, viscous advection-dominated accretion disk structure. 
The disk is tenuous because it accretes matter at well below the Eddington rate, and there- 
fore energy losses due to radiative processes are inefficient. General relativistic effects are 
taken into account in an approximate manner by utilizing the pseudo-Newtonian potential 
introduced by Paczyhski & Wiita (1980). The self-gravity of the disk is neglected. Under 
these assumptions, there are two conserved quantities in viscous ADAF disks, namely the 
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mass transport rate 

M = 4irrHpv, (1) 
and the angular momentum transport rate 

j = Mr 2 tt -g , (2) 

where r is the radial coordinate, H denotes the disk half-thickness, p is the mass density, 
v is the (positive) inflow velocity, f2 is the angular velocity, and Q denotes the torque. The 
energy transport rate is given by 

E = -gn + M^ + lvl + ^ + <^ , (3) 

where = r Q denotes the azimuthal velocity, P = (7 — 1) U is the pressure, U is the internal 
energy density, and $ represents the pseudo-Newtonian potential, defined by (Paczyhski & 
Wiita 1980) 

GM 

$ r = , (4) 

r-r s 

for a black hole with mass M and Schwarzschild radius r s = 2GM/c 2 . This potential 
successfully reproduces many aspects of the spacetime geometry around a nonrotating black 
hole. We assume that the adiabatic index, 7, maintains a constant value throughout the 
flow, and the transport rates M, J, and E are all defined to be positive for inflow. Under our 
assumptions, M and J are conserved throughout the disk, and therefore they represent the 
rates at which mass and angular momentum enter the black hole, respectively. In general, 
the energy transport rate E is also conserved, except at the location of an isothermal shock 
if one is present in the flow. 

The torque Q is associated with the gradient of the angular velocity Q through the 
standard relation (e.g., Frank, King, & Raine 2002) 

g(r) = -4nr 3 Hpu^ , (5) 
dr 

where v is the kinematic viscosity. In the present work, we adopt the Shakura-Sunyaev 
(1973) viscosity prescription, 

<r) = ^ , (6) 
where a is a constant, denotes the Keplerian angular velocity, defined by 

n U r) , °" = l -f , (7) 

r[r — r s ) z r dr 
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and a denotes the isothermal sound speed, defined by 

a 2 (r) = - p . (8) 

Combining equations ([!]), (J2J), ([5]), and (jSJ), we find that the gradient of the angular 
velocity can be written as 

dn_ vn K (e-£ Q ) 

J 9 9' \ ) 

dr ar z a z 

where £ = r 2 Q denotes the angular momentum per unit mass for the accreting gas at radius 
r, and £q = J/M is the angular momentum transport rate per unit mass. Since the torque 
Q vanishes at the event horizon (Becker & Le 2003), it follows from equation (j2J) that £q is 
also equal to the specific angular momentum of the material entering the black hole. The 
radial derivative of £ is given by 

d£ „ dQ , . 

* = 2r!i+r V (10) 

which can be combined with equation (JHD to obtain the differential equation 

d£ = 2£ v£ K (£-£o) , . 

7 9 9' \ ) 

ar r aa z r z 
where the Keplerian specific angular momentum £k is defined by 



r 3 / 2 VGM 

r) = r 2 n K (r) = . 12 

r — r 



s 



The disk half-thickness H is estimated using the standard vertical hydrostatic relation 

2 

nv 

H{r) = — . (13) 

In a steady state, the comoving radial acceleration rate in the frame of the accreting gas can 
be written as 

Dv dv ldP £l-£ 2 

= —v — = h — . (14) 

Dt dr p dr r 3 

By eliminating the torque Q between equations (j2J) and (EJ) and utilizing equation (IS]), we 
find that the energy transport rate per unit mass is given by 



E _ 1 2 1 £ 2 £ £ 7 
M 2 V 2 r 2 r 2 7 - 1 



e= — = ~ v ~ o^ + ^ + ^T^ + $ • ( 15 ) 



Since we are working within the ADAF framework, e will remain constant unless a shock 
is formed in the disk. In this scenario, the variation of the internal energy density U is 
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regulated by the adiabatic compression of the gas and the viscous dissipation in the disk. 
The comoving rate of change of U can therefore be written in the frame of the gas as (e.g., 
Becker & Le 2003), 

DU dU U dp pv (dl 2£\ 2 

-^ = -7-^ + ^1--- • (16) 



Dt dr p dr r 2 \dr r 

We can obtain the dynamical differential equation governing the spatial variation of the 
inflow velocity v by manipulating the conservation equations as follows. First we combine 
equations ([I]), f|T3|) . and (jSJ) to express the pressure as 

P = p±. (17) 

Using this relation to substitute for the pressure in the radial acceleration equation (TH|) . we 
find that the derivative of the sound speed a is given by 

da fa 2 \ v dv a d£^ 3a £ 2 — £%_ . . 

1 ) : — — — H h ' 1 ° 1 



dr \v 2 J a dr dr r ar 3 

Next we differentiate equation (IT51) and combine the result with equation ({TTI) to obtain 



dv £ 2 -l\ 2 7 a da e K (£- 1 v2 



(19) 



dr vr 3 7 — 1 v dr aa 2 r A 

By using equation (ITS!) to substitute for da/dr in equation (TT91 . we obtain the differential 
dynamical equation (e.g., Narayan et al. 1997) 

(± _ JlL\ l± = Jl_ (I _ + f7-l^ vW-io) 2 + t 2 -t 2 K (2Q) 

\a 2 y + ljv dr j + 1 \r £k dr J \y + 1 J aa 4 r 4 a?r 3 

In order to obtain the distributions of the flow variables v, a, and £, we must solve numerically 
the governing equations (ITT]) . (JT5|) . and (1201 . This requires consideration of the critical 
structure of the flow as discussed below. 



2.1. Critical Conditions 

A black hole will swallow matter either from its binary companion or from the winds 
of the surrounding stars. At the outer edge of the disk, the accreting matter has negligible 
radial velocity, although it ultimately crosses the event horizon with a local velocity equal 
to the speed of light. Hence accretion flows around black holes must be transonic. The 
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locations where the flow transitions from subsonic to supersonic are referred to as critical 
points. Depending on the parameters describing the properties of the gas at a large distance 
from the black hole, the flow may possess multiple critical points. When this occurs, it is 
interesting to investigate the possible existence of isothermal standing shock waves located 
between two of the critical points. The location of the shock, if one exists, is determined by 
applying the shock jump conditions. We first study the critical structure of the transonic 
flow by rewriting equation (120!) in the equivalent form 



dv _ N 
dr ~ ~D 1 

where the numerator and denominator functions N and D are defined by 



7 + 1 V 2r r — r c 



7 -i\ v 2 e K (e-e ) 2 , v(e 2 -£ 2 K ) 



and 



D 



7 + 1, 
v 2 27 



aa 4 r 4 



+ 



a 2 r 3 



a 2 7 + I 



(21) 



(22) 



(23) 



In order to obtain well-behaved global solutions, the numerator and denominator func- 
tions must vanish at exactly the same location in the flow. This condition represents the 
mathematical definition of the critical point. Setting N = D = 0, we obtain the critical 
conditions 



7 + 1 



3 

2rv 



1 



7-1 
7 + 1 



v 2 c e K (t c -to) 



+ 



a 2 .rl 







and 



2 7 







(24) 



(25) 



< 7 + 1 

where r c , v c , a c , and £ c denote the values of the radius, the velocity, the sound speed, and 
the specific angular momentum at the critical point, respectively. Any accretion solution 
must pass smoothly through the critical point, which is guaranteed if the critical conditions 
are satisfied. In this work, we assume approximate equipartition between gas and magnetic 
pressure following Narayan et al. (1997), and accordingly we set 7 = 1.5 in the subsequent 
analysis. 

The topological nature of a critical point is determined by the value of the velocity 
gradient dv/dr at that point (e.g., Das 2007; Becker & Le 2003). In general, we obtain two 
possible values for dv/dr, with one corresponding to accretion and the other to an outflow 
(wind). Any physically acceptable solution must pass through a saddle type or X-type critical 
point which is obtained when both of the derivatives are real and of opposite signs. If a shock 
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forms, then the accretion flow must also pass through another saddle type critical point in 
the post-shock region. When both of the values of dv/dr are complex, the critical point is 
referred to as O-type. This type of critical point is unphysical and therefore no acceptable 
flow solution can pass through such a point. 



2.2. Inner Boundary Conditions 

Solutions for the flow variables v(r), £(r), and a(r) are obtained by numerically integrat- 
ing equations ( TTTT) and ( 1201) . supplemented by the algebraic relation given by equation ( TT5l) . 
Becker & Le (2003) determined that in order to ensure the stability of the calculation, the 
integration must proceed in the outward direction, starting from a point close to the event 
horizon. This requires the availability of asymptotic boundary conditions that can be used 
to establish the starting values for the physical variables. Based on the fundamental general 
relativistic principle that the torque must vanish at the horizon (Weinberg 1972), Becker & 
Le (2003) derived the explicit asymptotic behaviors of v (r), £(r), and the "entropy function," 
K(r), for material subject to the Shakura-Sunyaev viscosity prescription, where 



K{r) 



7+1 



n 



K 



(26) 



The physical significance of K can be understood by combining equations ([I]), ffT31 . (jSJ), and 
( 126]) to show that 



#7-1 



oc 



(27) 



This relation establishes that K remains constant in regions of the flow unaffected by dissi- 
pation, where P oc p 7 . 



If the gas is in local thermodynamic equilibrium, then one can demonstrate that the 
value of K is related to the entropy per particle S by (Reif 1965) 



S = k In K + c 



(28) 



where cq is a constant that depends only the composition of the gas, but is independent of 
its state. Becker & Le (2003) established that the asymptotic behaviors of K(r) and £(r) 
close to the event horizon are given by 

1/2 



K r) = Kn 



16 a£ 2 
5 c 2 r^ 



- r W 2 



and 



2a. 



c 2 r c 



1/2 / K 2 \ £T 
^0 



2r 3 

s 



, 7+5 



S ' 



(29) 



(30) 
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respectively, where K$ represents the entropy value at the horizon. Note that the asymptotic 
radial dependence of K(r) is very weak, reflecting the unimportance of viscous dissipation 
near the horizon. Becker & Le (2003) also found that the asymptotic variation of the inflow 
velocity is given by 



v(r) = Vfi(r) 



1 + 



2eor 2 -^-( 7 + l)/(r) 
r 2 vl{r) - (7 - 1) f{r) 



1/2 



s ' 



(31) 



where eo denotes the value of the specific energy transport rate at the horizon, the free-fall 
velocity is given by 



v s {r 



r — r c 



(32) 



and the function / is defined by 



/(<•> - ^ 



r 



K 2 
-"■0 



2r 3 (r — r s ) 



7-1 
7+1 



(33) 



For given values of the parameters e , £0, and K , equations ([501 and (l3Tj) can be used to 
compute starting values for the integration of the variables £(r) and v(r) away from the event 
horizon. 



3. GLOBAL SOLUTIONS 

The structure of the global disk/shock model depends on the energy transport rate at 
the horizon eo, the angular momentum at the horizon £ Q , the entropy at the horizon K , 
the viscosity parameter a, and the ratio of specific heats 7. In the calculations performed 
here, we set 7 = 1.5, following Narayan et al. (1997). When a shock is present, we use the 
subscripts "-" and "+" to denote quantities measured just upstream and just downstream 
from the shock, respectively. In a shocked disk, the gas passes through one critical point 
in the pre-shock region at r = r° ut , and then through another in the post-shock region at 
r = r™ (Abramowicz & Chakrabarti 1990). Conversely, if the flow is smooth (shock- free), 
then the gas only passes through one critical point, located at radius r = r c . 

The process of determining the structure of a viscous ADAF disk containing an isother- 
mal shock begins with the selection of input values for e , £0, a, and 7. We must also select 
a provisional value for the entropy parameter K . Based on this information, we can utilize 
the asymptotic relations given by equations fl30l) and fl3Tl) to compute values for £{r) and 
v(r), respectively, at the starting radius for the outward integration, which is taken to be 
2.001 GM/c 2 . Next we numerically integrate equations ffTTj) . ffl5l) . and (TSUI) , beginning at the 
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starting radius and working outward towards the inner critical point. The value of Kq is var- 
ied until a solution is obtained that satisfies the critical conditions given by equations (jUj) 
and (|25p . which ensure that the flow passes smoothly through the inner critical point at 
radius r™. 

After the location of the inner critical point has been established, the integration must 
be continued into the upstream, subsonic region. Due to the mathematical structure of the 
dynamical equation (l2"Tj) . it is not possible to integrate across the critical point itself. We 
must therefore employ linear extrapolation, which requires the calculation of the derivative 
dv/dr at the critical point using FHopitaFs rule. After performing the linear extrapolation, 
the integration is continued in the subsonic region until we reach the shock radius, r*, which 
is assigned a provisional value initially. Since no energy is lost from the disk between the 
shock and the horizon, it follows that €q = e + , where e + is the energy transport rate on the 
downstream side of the shock. At the shock location, the velocity and the energy transport 
rate experience discontinuous transitions described by the isothermal shock jump conditions, 
which can be written as (e.g., Chakrabarti 1989a) 

V ^_=Ml, Ae = e + -e-= V -±^, (34) 

where Ai + = v + /a + denotes the post-shock Mach number. This relation facilitates the 
calculation of the pre-shock (incident) energy transport rate, e_. Note that if no shock is 
present, then e is continuous throughout the entire flow, and therefore Ae = and e = e_. 
By definition, the flow maintains a uniform temperature across the isothermal shock, and 
consequently a + = a_. The shock radius r* is varied until we obtain a solution that passes 
smoothly through the outer critical point, located in the upstream region at radius r° ut . The 
solution process is completed when the critical conditions are satisfied at the outer critical 
point. For given input values of eo, io, a, and 7, in the end we obtain unique values for K , 
r*, r™, and r° ut , in addition to the associated global dynamical solution for the disk/shock 
structure. 

A typical set of shocked and smooth accretion solutions obtained for the same value 
of the incident energy transport rate e_ is presented in Figure la. In these calculations, 
the Shakura-Sunyaev viscosity parameter was specified by setting a = 0.1, and the specific 
heats ratio is given by 7 = 1.5. By applying the multi-step iterative method described 
above, we obtain for the shocked flow the parameters e_ = 0.001856, e = e + = —0.003, 
£0 = 2.55, and K = 0.00689. In the shocked case, the incident flow first becomes supersonic 
upon passing through the outer critical point located at r° ut = 73.51. The gas subsequently 
passes through an isothermal shock at r* = 31.2, where the velocity makes a discontinuous 
transition to become subsonic again. Downstream from the shock, the flow accelerates 
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10 1 10 2 10 3 

r 

Fig. 1. — (a) Inflow velocity v (solid lines) and isothermal sound speed a multiplied by 
a/27/ (7 + 1) (dashed lines) plotted as functions of the radius r in units of GM/c 2 for a 
typical global accretion solution with 7 = 1.5. The thick and thin lines represent the shocked 
and smooth solutions, respectively, (b) Energy transport rate e plotted as a function of radius 
r. Note the jump to the negative post-shock value due to the release of energy through the 
upper and lower surfaces of the accretion disk at the isothermal shock location. See the text 
for details. 
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Fig. 2. — Specific angular momentum t plotted as a function of the radius r for (a) a = 0.01, 
(6) a = 0.1, and (c) a = 0.2. Here we assume that 7 = 1.5. The solid, dotted, and dashed 
lines correspond to the shocked, smooth, and Keplerian profiles, respectively. See the text 
for the complete sets of parameters used in each panel. 
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towards the horizon, and regains its supersonic character after crossing the inner critical 
point at r™ = 6.798. The radius at the outer edge of the hot region, r edge = 3555, is 
computed by setting £ = £^, so that the angular momentum matches the local Keplerian 
value. The radial profile of the specific angular momentum £ is plotted in Figure 2b. Although 
the incident energy transport rate e_ is negative in our calculations, it should be emphasized 
that the total energy inflow rate, including the rest mass energy, is positive as required. 

For the smooth- flow solution depicted in Figure la, we used the parameters e = e + = 
e_ = 0.001856, £q = 2.368, and K = 0.008565. In this case, the value of £q was selected 
so that £ = £k at the same outer edge radius, r c d gc = 3555, obtained in the shocked disk 
(see Fig. 26). The smooth flow solution passes through a single critical point, located at 
r c = 8.097. It is interesting to compare the plots of the sound speed obtained in the smooth 
and shocked cases. The loss of energy at the shock significantly reduces the sound speed 
(and therefore the temperature) in the post-shock region relative to the smooth solution. In 
Figure lb, we plot the variation of the energy transport rate for the shock solution. The 
excess energy released through the upper and lower surfaces of the disk at the shock location 
may be sufficient to power the formation of a relativistic jet (Le & Becker 2005; Becker, Das, 
& Le 2008). 

The viscosity parameter a plays a central role in determining the efficiency of the angular 
momentum transport in the disk. We explore the effect of varying this parameter in Figure 2, 
where we plot the results obtained for the radial profile of the specific angular momentum 
£ in both shocked and smooth disks for three different values of a. All of the solutions 
are sub-Keplerian (£ < £■&) for r < r c d gc . In the shocked case, the increase in the angular 
momentum in the post-shock region produces the "centrifugal barrier" required for shock 
formation, while no such barrier exists in the smooth case. In each of the panels, the value 
of £o was determined by satisfying the inner critical conditions, and the incident energy 
transport rate e_ has the same value in the shocked and smooth scenarios. The parameter 
values adopted in panel (a) are a = 0.01, 7 = 1.5, e + = —1.0 x 10 -3 , e_ = 1.9072 x 10~ 3 , 
£ = 3.0, K = 2.7509 x 10" 3 in the shocked case, and e+ = e_ = 1.9072 x 10" 3 , £ = 2.845, 
K = 3.8588 x 10~ 3 in the smooth (shock-free) case. In panel (6), we adopt the same 
parameters used in Figure 1, namely a = 0.1, 7 = 1.5, e + = —3.0 x 10~ 3 , e_ = 1.856 x 10~ 3 , 
£ = 2.55, K Q = 6.89 x 10~ 3 in the shocked case, and e+ = e_ = 1.856 x 10" 3 , £ = 2.368, 
K = 8.565 x 10~ 3 in the smooth case. The values adopted in panel (c) are a = 0.2, 7 = 1.5, 
e+ = -1.0 x 10~ 3 , e_ = 1.6483 x 10~ 3 , £ = 2.1, K = 1.1883 x 10" 2 in the shocked case, 
and e + = e_ = 1.6483 x 10" 3 , £ = 1.966, K = 1.25 x 10~ 2 in the smooth case. Note the 
increasing prominence of the centrifugal barrier as a increases. 
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4. SHOCK PROPERTIES 

Since this is the first work to extensively consider the consequences of a standing shock 
for the structure of a viscous accretion disk, it is interesting to examine how the shock 
properties depend on the various flow parameters. We carry out the analysis for flows char- 
acterized by the fundamental parameters (e , £o, a), where e denotes the energy transport 
rate at the horizon, which is the same as the post-shock energy transport rate e+. In each 
case, we keep one of the three parameters (eo, £o, &) fixed and vary the other two to study 
the resulting shock properties. 

In Figure 3 we present the results obtained for the shock radius, r*, the compression 
ratio, R = V-/v + , and the energy jump, Ae = e + — e_, as functions of the angular momentum 
transport rate £ and the energy transport rate at the horizon e . In these calculations, 
the viscosity and specific heat parameters are held constant, with a = 0.1 and 7 = 1.5, 
respectively. Figure 3 a clearly demonstrates that accretion flows with a very wide range of 
input parameters may possess isothermal standing shocks. Note that for a given post-shock 
energy parameter eo, the shock location recedes away from the black hole with increasing £ 
because of the enhancement of the centrifugal barrier. On the other hand, the shock location 
shifts inward as eo decreases for a given value of £0, due to the increased ram pressure of the 
flow. For a given value of e , there is a restricted range of £ Q values within which a shock can 
exist. In general, we find that the shock radius r* always exceeds ~ 10 GM/c 2 . In Figure 36, 
we plot the variation of the shock compression ratio R as a function of the flow parameters 
e and £q. We observe that for fixed eo, the compression ratio increases as the shock moves 
away from the black hole. In Figure 3c, we plot the corresponding variation of the energy 
jump Ae. Note that Ae becomes more negative as eo decreases for fixed £q, which reflects 
the fact that a larger amount of gravitational potential energy is released into the outflow 
when the shock wave forms closer to the black hole. 

Next we examine the dependence of the shock properties on the viscosity parameter a. 
In Figure 4 a, the isothermal shock location r* is plotted as a function of a for various values 
of the post-shock (horizon) energy transport rate e for the case with 7 = 1.5. Here, we fix 
the angular momentum transport rate by setting £ = 2.62. In all cases, we observe that 
shocks can form over a wide range of flow parameters. For a given a, the shock tends to form 
farther away from the black hole as e is increased, which is consistent with the behavior 
noted in Figure 3a. We also observe that for fixed e , the shock radius r* increases with 
a. This is because higher values of a tend to reduce the angular momentum, which forces 
the centrifugal barrier to operate at a larger radius in order to be effective. Note that for 
a given set of flow parameters (eo, £0), the viscosity parameter a has both upper and lower 
limits for shock formation. The lower limit reflects the fact that fully inviscid flows cannot 
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Fig. 3. — Variation of (a) the isothermal shock radius r*, (b) the compression ratio R = 
v-/v+, and (c) the energy jump Ae as functions of £ . The solid, dotted, dashed, dot- 
dashed, and short-long-dashed curves correspond to e = —0.006, —0.005, —0.004, —0.003, 
and —0.002, respectively. Here, we set a = 0.1 and 7 = 1.5. 




Fig. 4. — Variation of (a) the isothermal shock radius r*, (b) the compression ratio R = 
v-/v+, and (c) the energy jump Ae as functions of a. The solid, dotted, and dashed lines 
denote e = —0.006, —0.005, and —0.004, respectively. Here, we set £ = 2.62 and 7 = 1.5. 
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form shocks for arbitrary combinations of e and £ (e.g., Le & Becker 2005). Conversely, 
the upper limit on a is due to the efficient dissipation of angular momentum that occurs 
as a increases. This effect prevents the development of the centrifugal barrier required for 
shock formation if a is too large. 

The dependence of the shock properties on post-shock energy e is illustrated in Figure 5 
for three values of the angular momentum transport rate £q. The viscosity and specific heat 
parameters are held constant with the values a = 0.1 and 7 = 1.5, respectively. In Figure 5a 
we plot the shock location r* as a function of eo- Since the standing shock is centrifugally 
supported, the shock radius r* increases when £q is increased for a fixed value of eo- On the 
other hand, if we hold r* constant, then we observe an anticorrelation between £0 and eo- In 
Figure 56 and Figure 5 c, respectively, we depict the associated variations of the compression 
ratio R and the energy jump Ae as functions of e for the same set of flow parameters used 
in Figure 5 a. The compression ratio increases when either £ or eo is increased. We note 
that Ae increases almost linearly with increasing eo for fixed £q. The results indicate that 
flows with larger values of £q tend to discharge less energy at the shock location. 

The hot, sub-Keplerian ADAF region must connect with a cool Keplerian outer disk at 
some large radius. We estimate the outer edge radius of the ADAF zone, r cdgc , by setting the 
angular momentum of the accretion flow equal to the local Keplerian value. In Figure 6 we 
plot the variation of the outer edge radius as a function of £q and eo for the case with a — 0.1 
and 7 = 1.5. We observe that for a fixed value of e , the value of r e d ge rapidly increases with 
increasing £ , indicating that the cool Keplerian disk recedes from the horizon. The outer 
edge radius also increases with increasing eo when £q is held fixed. 

The shocked, viscous, advection-dominated disks discussed here represent a new type 
of astrophysical phenomenon, and therefore it is important to investigate the range of flow 
parameters within which shocks can form. In Figure 7 we analyze this question for the 
case with 7 = 1.5 by varying the angular momentum transport rate £q and the upstream 
(incident) energy transport rate e_ while holding the viscosity parameter a fixed at five 
different values. The solid curve depicts the region within which shocked solutions are 
possible for fully inviscid flow (a — 0). This result agrees very well with Figure 2 from Le & 
Becker (2005). We note that the allowed region for the shocked solution shifts toward lower 
angular momentum as a increases. This shift occurs because the specific angular momentum 
of the accreting material at the horizon, £0, decreases in response to the enhanced viscosity 
associated with the increase in a. The shrinkage of the region for the shocked solution 
observed as a increases suggests that shocked disk solutions are impossible above a critical 
value of a, denoted by a C rit- For the model considered here, with 7 = 1.5, we find that 
a crit « 0.27. 
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Fig. 5. — Variation of (a) the isothermal shock radius r*, (b) the compression ratio R = 
v-/v+, and (c) the energy jump Ae as functions of eo- The solid, dotted, and dashed curves 
represent £ = 2.62, 2.67, and 2.72, respectively. Here, we set a = 0.1 and 7 = 1.5. 
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Fig. 6. — Plot of the outer edge radius r e d ge as a function of £q- This is the radius where 
the angular momentum £ equals the Keplerian value £k (eq. [12] )• The solid, dotted, and 
dashed lines correspond to e = 0.0, —0.004, and —0.008, respectively. Here, we set a = 0.1 
and 7 = 1.5. See the text for details. 
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Fig. 7. — Representation of the parameter space spanned by the specific angular momentum 
at the horizon £ and the incident energy transport rate e_. The wedge-shaped regions denote 
the areas of the parameter space within which shocks can form for the indicated values of 
a. See the discussion in the text. 
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5. ASTROPHYSICAL APPLICATIONS 

It is interesting to ask whether the self-consistent disk/shock model we have described 
can account for the energetics of specific jet sources, such as M87 and Sgr A*. When observa- 
tional values of the black hole mass M, the accretion rate M, and the jet kinetic luminosity 
Lj et are available for a particular source, we can utilize this information to constrain the 
model parameters. The process begins by selecting values for the fundamental free parame- 
ters eo, a, and 7. After selecting provisional values for C.q and Kq, we utilize the asymptotic 
relations given by equations (1301 and (I3TI) to compute initial values for the functions £{r) 
and v(r) at the starting radius 2.001 GM/c 2 . We then numerically integrate equations fllip . 
f ll5p . and (1201) in the outward direction, and vary the value of K Q until the critical conditions 
given by equations (1241) and (1251) are satisfied at the inner critical point. 

The energy transport rate e drops from its upstream value e_ to the downstream value 
e + = eo at the shock location due to the escape of energy from the disk into the outflow (jet). 
We can therefore determine the shock radius r* by requiring that the value of Ae computed 
using equation (1341) agrees with that calculated based on the observed jet kinetic luminosity, 
Lj et , using the relation 

L jet = -MAe , (35) 

where the negative sign appears because Ae < 0. Once the shock radius is established, the 
integration is continued on the upstream side of the shock toward the outer critical point. 
If the flow does not pass smoothly through the outer critical point, then the value of £q 
is modified and the procedure is repeated starting from the horizon. The computation of 
the self-consistent disk/shock dynamical model is complete when the critical conditions are 
satisfied at the outer critical point. For given input values of e , a, and 7, in the end we 
obtain the complete global dynamical solution along with unique values for £ , K , r™, 
and r° ut . We provide two specific examples below. 

5.1. M87 

The disk/shock/outflow model we have developed was previously utilized by Becker, 
Das, & Le (2008) to analyze the structure of the accretion flow in M87. The results obtained 
for that source are presented in more detail here. Our application to M87 is based on 
the observational values M = 1.34 x lO^ 1 M yr _1 (Reynolds et al. 1996), Lj Ct = 5.5 x 
10 43 erg s" 1 (Reynolds et al. 1996; Bicknell & Begelman 1996; Owen, Eilek, & Kassim 2000), 
and M = 3.0 x 10 9 M Q (e.g., Ford et al. 1994). We set the viscosity parameter a = 0.1 in 
order to demonstrate that shocks can exist in ADAF disks even in the presence of substantial 
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viscosity, although our model can accommodate any value for a. Based on the assumption 
of approximate energy equipartition between the magnetic field and the gas internal energy, 
we adopt the value 7 = 1.5 for the adiabatic index (Narayan et al. 1997). We utilize natural 
gravitational units in our numerical examples, with GM = c = 1 and r s = 2, except as 
noted. The remaining parameters for the shocked-disk model can be computed based on the 
observations of M87, from which we obtain e_ = 0.001516, e+ = e = -0.005746, £ = 2.6257, 
K = 0.00608, r* = 26.329, = 6.462, and r° ut = 96.798. For the pre-shock and post- 
shock velocities and Mach numbers we obtain w„ = 0.138, u + = 0.068, M._ = 1.427, and 
M. + = 0.701, respectively, and for the disk half-thickness at the shock location we obtain 
H, = 12.10. 

The results obtained for the inflow speed u and the isothermal sound speed a are plotted 
in Figure 8a, and the associated solution for the specific angular momentum £ is depicted 
along with the Keplerian profile £k (eq. [12]) in Figure 8b. Results are presented for both 
shocked and smooth (shock-free) disks. In the shocked disk, the source is located at the 
shock radius, r* = 26.329, and we observe that the values of £ and £k merge at r e d ge = 4658, 
which represents the outer edge of the ADAF region. The smooth solution is computed by 
integrating the conservation equations using the same value for the incident energy transport 
rate e_ as adopted in the shocked-disk model, but with no shock included. The value of £ is 
then varied until we obtain £ = £k at the same outer radius, r e d ge = 4658, as in the shocked 
disk. The resulting parameter values for the smooth solution are e_ = e + = 0.001516, 
£q = 2.3988, K = 0.0084, and r c = 7.572. In order to facilitate comparison with the shocked 
model, in the smooth case the source is assumed to be at the same radius r = 26.329. 

Our solutions for the specific angular momentum £ are significantly sub-Keplerian, in 
agreement with the results obtained by Narayan et al. (1997), who utilized the same set of 
ADAF conservation equations employed here. The sub-Keplerian nature of the flow stems 
from the relatively high value for the viscosity parameter adopted here, a = 0.1. Note that 
the sound speed (and hence the temperature) is significantly lower in the shocked case as 
compared with the smooth disk due to the release of energy into the outflow at the shock 
location. 

It is important to examine the energy balance in the shocked model to ensure that 
the rate of energy loss at the shock location is equal to the observational value for the 
kinetic luminosity in M87. We can compute the theoretical value for the jet luminosity, 
Lj et , by substituting our values for e + and e_ into equation (|35|) . The result obtained is 
Lj e t = 5.5 x 10 43 erg s _1 , which agrees with the observed jet luminosity for M87. This 
confirms that our viscous shocked disk model is able to account for the energetics of the M87 
outflow. 
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Fig. 8. — Dynamical profiles for M87. (a) Inflow velocity v (solid lines) and isothermal 
sound speed a multiplied by a/^t/vT + -0 (dashed lines) plotted as functions of the radius 
r in units of GM/c 2 . The thick and thin lines denote the shocked and smooth solutions, 
respectively, (b) Specific angular momentum I for the shocked (thick line) and smooth (thin 
line) solutions plotted as functions of r in units of GM/c 2 along with the Keplerian angular 
momentum £k (dot-dashed line). 
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5.2. SgrA* 

For SgrA*, we adopt the values M = 2.6 x 10 6 M Q (Schodel et al. 2002) and M = 
8.96 x 10~ 7 M yr" 1 (Yuan, Markoff, & Falcke 2002; Quataert 2003), and we set 7 = 1.5 
and a = 0.1. The kinetic luminosity of the jet is estimated using L )et = 5 x 10 38 ergs s _1 
(Falcke & Biermann 1999), although this value is rather uncertain (e.g., Yuan 2000; Yuan et 
al. 2002). The results for the remaining model parameters implied by the observations are 
e_ = 0.00134884, e+ = e = -0.0085, £ = 2.6728, K = 0.005448, r* = 19.917, rf = 6.380, 
and r° ut = 112.384. For the pre-shock and post-shock velocities and Mach numbers we 
obtain m_ = 0.159, u + = 0.0748, M- = 1.4578, and M+ = 0.6860, respectively, and for the 
disk half-thickness at the shock location we obtain = 8.72. 

In Figure 9 a we plot the results obtained for the inflow speed u and the isothermal 
sound speed a in both the shocked and smooth disk models for SgrA*. The profiles of £ and 
£k are plotted in Figure 9b. In the shocked case, the source is located at r* = 19.917, and 
we find that £ = £k at radius r edge = 5432, which is the outer edge of the ADAF region. The 
incident energy transport rate e_ = 0.00134884 is the same for both the smooth and shocked 
models, and the value of £q in the smooth case is determined by requiring that £ merge 
with the Keplerian profile at the same radius r e d ge = 5432 as in the shocked case. Based 
on this approach, the parameters obtained for the smooth model are e_ = e + = 0.00134884, 
£q = 2.416, K = 0.008293, and r c = 7.511. We assume for consistency that the source is 
located at the same radius r = 19.917 in both the smooth and shocked cases. 

The results for the specific angular momentum £ are significantly sub-Keplerian, as 
expected based on the relatively high value utilized for the viscosity parameter, a = 0.1 
(cf. Fig. 86). The sound speed is lower when a shock is present, due to the escape of 
energy at the shock location, as noted previously in Figure 8a. The result for L- ]et obtained 
by substituting our values for e + and e_ into equation (|35|) is equal to the observed jet 
luminosity, 5 x 10 38 erg s -1 , and this establishes that our shocked disk model successfully 
accounts for the observed energetics of the accretion flow and the outflow in SgrA*. 

6. CONCLUSION 

In this paper we have presented the first systematic study of the structure of stationary, 
advection-dominated, shocked viscous accretion disks around black holes. Our primary focus 
here has been on the global properties of accretion flows containing isothermal shock waves, 
formed as a result of the presence of a centrifugal barrier located near the event horizon. 
Our model is based on the same set of standard ADAF conservation equations considered by 
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Narayan et al. (1997). The conservation equations are supplemented by the inner boundary 
conditions developed by Becker & Le (2003) which allow the calculation of starting values 
for the physical variables close to the horizon. Once the starting values are established, a 
stable outward integration method is employed to obtain the radial profiles of the various 
flow quantities such as velocity, sound speed, and angular velocity. The fundamental model 
parameters are the energy inflow rate at the horizon, eo, the specific angular momentum 
of the material at the horizon, £ , the entropy parameter at the horizon, K Q , the Shakura- 
Sunyaev viscosity parameter, a, and the ratio of specific heats, 7. 

For given observational values of the black hole mass M, the accretion rate M, and the 
jet kinetic luminosity Lj et , we are free to select values for a, 7, and eo- Following Narayan et 
al. (1997), we have set 7 = 1.5 in the calculations presented here. The remaining parameters 
£ and K , and the shock radius r*, are subsequently determined by applying the asymptotic 
inner boundary conditions given by equations (1301) and (l3~Tj) . the critical conditions given by 
equations IH)) and (I25I) . and the shock jump conditions given by equations (1341) and (|35|) . 
The dynamical structure of the disk/shock system is obtained by numerically integrating 
equations ( TTT1) . (fl5l) . and (1201 in the outward direction beginning at a starting radius located 
close to the event horizon. 

The formation of a shock requires that the flow possess at least two critical points. 
The existence of multiple critical points in both adiabatic and polytropic accretion flows 
has been reported by many authors (e.g., Abramowicz & Zurek 1981; Chakrabarti 1989a, 
1989b, 1990, 1996). Flows with multiple critical points are confined to a restricted region of 
the parameter space, as indicated in Figure 7. Accretion flows within this region can form 
isothermal shocks. In their earlier study, Narayan et al. (1997) focused solely on solutions 
that pass only through the inner critical point, and consequently they did not treat the 
shocked disk solutions. Flows with discontinuities must first pass through an outer critical 
point, and subsequently through an isothermal shock and an inner critical point in order to 
represent physically acceptable accretion solutions. That is the scenario we have investigated 
here. 

The possible presence of shock waves in accretion disks and their observational impli- 
cations have been examined by a variety of previous authors. Abramowicz & Chakrabarti 
(1990) presented one of the earliest studies of the properties of standing shock waves in accre- 
tion disks. Subsequently, the particular importance of isothermal shock waves was realized 
by Le & Becker (2004, 2005, 2007), who pointed out that efficient particle acceleration in the 
vicinity of a standing isothermal shock can help to explain the formation of the relativistic 
outflows of matter commonly observed around radio-loud, underfed black holes. However, 
all of these previous studies of isothermal shocks in accretion disks neglected the important 
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role of viscosity in determining the disk structure. Very recently, Becker, Das, & Le (2008) 
demonstrated that first-order Fermi acceleration inside a viscous, shocked ADAF can ac- 
celerate particles more efficiently than in a smooth disk. Hence the presence of a standing 
shock wave can create a potentially favorable environment for the production of relativis- 
ts outflows. In the present paper, we have presented a detailed study of the dynamics of 
viscous, shocked ADAF disks that are governed by the standard description of the angular 
momentum transport based on the Shakura-Sunyaev prescription for the kinematic viscosity. 
In particular, we confirm that the energetics of the viscous disk/shock model can explain 
the observational properties of both the accretion flow and the relativistic outflows in M87 
and Sgr A*. 

In Figure 7 we explored the parameter space within which ADAF disks with isothermal 
shocks can exist as a function of the specific angular momentum at the horizon £ Q and the 
incident energy transport rate e_ for various (constant) values of the viscosity parameter 
a. The shock region shifts toward lower values of £o as a increases due to the enhanced 
outward diffusion of angular momentum. We find that for the canonical value 7 = 1.5, 
shocked disk solutions can exist provided a < 0.27. This establishes for the first time that 
highly viscous ADAF disks can contain standing shock waves. The stability of such a shock 
is an open question, but Nagakura & Yamada (2008) and Okuda et al. (2007) recently 
found that shocks in inviscid disks are persistent, although the shock radius may oscillate. 
These results at least suggest the possibility that shocks in viscous disk may be stable as 
well, although this needs to be investigated in future work. Since shocked solutions are 
expected to possess higher entropy than smooth solutions, we argue based on the second 
law of thermodynamics that when shock solutions are permitted dynamically, they should 
be expected to form (Becker & Kazanas 2001). 

Although Gu & Lu (2001, 2004) obtained solutions for ADAF disks containing Rankine- 
Hugoniot shocks, the results developed here represent the first dynamical solutions for ADAF 
disks with a significant level of viscosity (i.e., a = 0.1) containing isothermal shocks. The 
Rankine-Hugoniot shocks studied by Gu & Lu (2001, 2004) have conserved energy transport 
rates, and therefore they cannot directly produce relativistic outflows (jets). Conversely, the 
isothermal shocks we study here have a discontinuous energy transport rate and consequently 
they are able to power the observed outflows. In our analysis of the disk structure, we have 
assumed that the accretion rate M remains constant throughout the disk, and therefore we 
have neglected the effect of the outflow mass loss on the dynamics of the disk, despite the 
fact that an outflow emanates from the shock location in our model. This assumption is 
probably reasonable since estimates show that the rate of mass loss from the shock into the 
wind is typically ~ 10^ 3 M if the jet is hadronic (Le & Becker 2005); this ratio is further 
reduced by the factor m e /m p if the jet is leptonic. However, there are two other potential 
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forms of mass loss from the disk that must also be considered. The first is the possibility 
of hydrodynamically driven outflows of the background thermal gas, and the second is the 
possibility of the escape of a significant number of relativistic particles from the inner and/or 
outer regions of the disk, rather than from the shock itself. We consider each of these issues 
below. 

Narayan et al. (1997) and Blandford & Begelman (1999) have noted that ADAF disks 
typically possess a positive Bernoulli parameter, 

„, , v 2 I 2 7 a 2 GM . . 

2 2r z 7 — 1 r — r s 

suggesting that sufficient energy is available to power hydrodynamical outflows of the thermal 
gas. It is interesting to reexamine this issue in the context of the shocked, viscous accretion 
disks considered here. In Figures 10a and 106, respectively, we plot the Bernoulli parameter 
as a function of radius based on the smooth and shocked models for M87 and Sgr A* developed 
in § 5. Both the smooth and shocked solutions display a slightly positive Bernoulli parameter 
in the inner and outer regions, suggesting that the gas is unbound there. When a shock is 
present, the gas is bound in the post-shock region, but it becomes unbound again before 
crossing the event horizon. 

The Bernoulli parameter is related to the energy transport rate E and the viscous torque 
Q via (see eq. [3]) 

E = —Q Q + MB . (37) 

Since E is conserved in the inner region, it follows that the positive bump in the Bernoulli 
parameter between the shock and the horizon reflects an increase in the torque-driven out- 
ward energy transport rate, Q Q. This effect is less pronounced in the shocked solution 
because the escaping particles carry away a significant fraction of the binding energy. The 
positivity of the Bernoulli parameter in the inner and outer regions suggests the possibility 
of hydrodynamical outflows emanating from the non-shock regions of the flow. However, 
any such outflows would be quite weak, with terminal Lorentz factors T = 5 + 1 ~ 1.01, 
whereas the terminal Lorentz factor for the outflow from the shock is T ~ 6 (Becker et al. 
2008; 2009). More fundamentally, it must be emphasized that despite the positivity of the 
Bernoulli parameter, there is no hydrodynamical mechanism in the model considered here 
capable of generating an outflow, since the ADAF disks we focus on are in vertical hydro- 
static equilibrium, and therefore the net vertical force vanishes (cf. eq. [IS])- Hence some 
additional mechanism, outside the standard ADAF framework, must be operative in order 
to channel the energy into the outflows. In our view, the outflows are not hydrodynamical in 
nature, but are instead powered by the first-order Fermi acceleration of relativistic particles 
in the vicinity of the discontinuous shock. 
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In addition to the possibility of hydrodynamical mass loss, some mass may also be lost 
from the disk in the form of relativistic particles that escape from the inner and out regions of 
the disk, rather than from the shock itself. The escape of these additional relativistic particles 
could be driven by the distributed first-order Fermi acceleration occurring throughout the 
disk as a consequence of the overall compression of the flow. The significance of this effect 
can only be understood by performing a detailed study of the relativistic particle transport 
occurring throughout the disk, which is beyond the scope of the present paper. We are 
currently developing such a model and will report the results in a separate paper. Preliminary 
calculations suggest that the spatial distribution of the escaping particles is strongly peaked 
at the shock location, with about half of the particles escaping from the shock itself, and the 
remainder escaping from the surrounding region (Becker et al. 2009). Hence, we conclude 
that the loss of mass from the disk is strongly concentrated in the vicinity of the shock, and 
that it is dominated by the escape of relativistic particles, rather than by hydrodynamical 
outflows of the background (thermal) gas. 

The authors would like to thank the anonymous referee for providing several useful com- 
ments that significantly strengthened the paper. SD acknowledges support via a postdoctoral 
fellowship from the Korea Astronomy and Space Science Institute (KASI). 
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Fig. 9. — Same as Fig. 8, except model results are for Sgr A*. 




Fig. 10. — Bernoulli parameter (eq. [36]) plotted as a function of the radius r in units of 
GM/c 2 for (a) M87 and (b) Sgr A*. The thick and thin lines denote the shocked and smooth 
solutions, respectively. 
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